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The random forced Navier-Stokes equation can be obtained as a variational problem of a proper 
action. By virtue of incompressibility, the integration over transverse components of the fields allows 
to cast the action in the form of a large deviation functional. Since the hydrodynamic operator is 
nonlinear, the functional integral yielding the statistics of fluctuations can be practically computed 
by linearizing around a physical solution of the hydrodynamic equation. We show that this procedure 
yields the dimensional scaling predicted by K41 theory at the lowest perturbative order, where the 
perturbation parameter is the inverse Reynolds number. Moreover, an explicit expression of the 
prefactor of the scaling law is obtained. 

I. INTRODUCTION 

A field theoretic approach to the study of the random stirred Navier-Stokes equation (rsNSE) can be traced back 
to the seminal paper by Martin, Siggia and Rose [1]. This was the starting point for the apphcation of many field- 
theoretic strategies, e.g. diagramatic expansions, renormalization group methods [2] (for recent developments and 
apphcations the reader can be addressed to [3] ) , instanton-based approaches (for applications of instantonic methods 
in turbulence see, e.g., [4-7] and references therein) and combinations of them [8]. The many technical difficulties 
encountered in developing these approaches avoided to gather conclusive achievements. 

In this paper we show that one step forward along this field-theoretic approach allows one to cast the action associated 
with the rsNSE into the form of a large deviation functional. Recently, large-deviation theory has scored sensible 
success in describing fiuctuations in stationary non-equilibrium regimes of various microscopic models [9]. This 
approach is mainly based on the extension of the time-reversal conjugacy property introduced by Onsager and Machlup 
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[10] to stationary non-equilibrium states. In practice, thermal fluctuations in irreversible stationary processes can be 
traced back to a proper hydrodynamic description derived from the microscopic evolution rules. The general form of 
the action functional is 



where p{t, x) represents in general a vector of thermodynamic variables depending on time t and space variables x. 
The symbol (• , •) denotes the integration over space variables. W is a hydrodynamic evolution operator acting on 
p: it vanishes when p is equal to the stationary solution p, which is assumed to be unique. The positive kernel K{p) 
represents the stochasticity of the system at macroscopic level. According to the large deviation-theory, the entropy 
S of a, stationary non-equilibrium state is related to the action functional / as follows: 



where the minimum is taken over all trajectories connecting p to p. 

For our purposes it is enough to consider that the action functional / provides a natural measure for statistical 
fluctuations in non-equilibrium stationary states, so that, formally, any statistical inference can be obtained from I. 
Indeed, from the very beginning we have to deal with a hydrodynamic formulation, namely the rsNSE: in the next 
Section we will argue that an action functional of the form (1) can be obtained by field-theoretic analytic calculations. 

In particular, explicit integration over all longitudinal components of the velocity field and over the associated 
auxiliary fields can be performed. This allows to obtain a hydrodynamic evolution operator W which depends only 
on the transverse components of the velocity field Vj.{t,x) {a = 1,2,3). Moreover, the positive kernel K amounts to 
the inverse correlation function of the stochastic source. This formulation allows to overcome some of the technical 
difficulties characterizing standard perturbative methods and diagramatic expansions. 

On the other hand, we have to face with new difliculties. The hydrodynamic operator appearing in the large 
deviation fimctional is nonlinear, so that functional integration is unfeasible. One has to identify a solution v^{t,x) 
of the associated hydrodynamic equation and linearize the hydrodynamic operator around such a solution. Then, 
functional integration can be performed explicitly on the "fluctuation" field. In order to be well defined, this ap- 
proximate procedure would demand the uniqueness of the solution of the nonlinear hydrodynamic equation. For this 
reason we have restricted our choice to a class of space-time functions which are also solutions of the linear problem. 
Among them, there is only one function which satisfies physically relevant boundary conditions (see Section III). 
Statistical fluctuations have been estimated with respect to this solution, which has also the advantage of reducing 
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S{p) =inf/[_<^^o](p) 
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the dependence of the generating functional on the pressure field to a trivial constraint. In practice, we construct a 
perturbative saddle-point approach based on a linearization procedure of the velocity field Vj<{t,x) around Vj,{t,x). 
As a consequence of the nonlinear character of the original problem, the fluctuation fleld x) = Vj:{t, x) — v^{t, x) 
is found to obey a linearized hydrodynamic problem with coefficients depending on space and time through Vj,{t,x). 
It is worth stressing that even the solution of the linearized problem is nontrivial and it is found to depend naturally 
on a perturbative parameter TZ~^, the inverse of the Reynolds number. We exploit this property by constructing a 
further perturbation procedure to obtain an explicit expression for u^{t,x) at different orders in TZ~^. These points 
are discussed in Section IV. 

Since our main purpose here is the estimation of the structure function (see Section V) as an average over the 
non-equilibrium measure induced by the action I, we have to assume that the perturbative expansion applies in a 
wide range of values of TZ. In particular, we guess that it holds also for moderately large TZ, since a statistical average 
of any observable cannot be valid for too large values of TZ, i.e. in a regime of fully developed turbulence. We will 
argue that statistical estimates can be consistently obtained for values of TZ which extend up to the region of stability 
of the solution Vj:{t, x). Beyond this region we have no practical way of controlling the convergence of the linearization 
procedure. It is worth stressing that we obtain an analytic expression of the structure function: the so-called K41 
scaling law [11] is recovered on a spatial scale, whose nontrivial dependence on TZ is explicitly indicated. 

At the present stage, we are not able to say at which extent our results on the dimensional scaling are dependent 
on the particular choice we did for the solution around which wc studied the fluctuations. Further investigations are 
needed to clarify this important point, which probably require the combination of analytical and numerical techniques. 



II. THE MODEL 

We consider the Navier-Stokes equation for the velocity vector-field components v°' {t, x) describing a divergence- free 
homogeneous isotropic flow: 

=0. (4) 

dx" ^ ^ ^ 

Here, P is the pressure and the field represents a source/sink of momentum necessary to maintain velocity 
fluctuations. Customarily [12], we assume /" to be a white-in-time zero-mean Gaussian random force with covariance 

{r{t,x)f{t',x'))=F'^''{x-x')6{t-t') . (5) 



Due to constraint (4), the field v°'{t,x) depends only on the transverse degrees of freedom of f°^{t,x). Without 
prejudice of generality we can also assume divergence-free forcing, yielding the additional relation 



A standard choice for F°'^ is 



(6) 



(7) 



where Z^o is the power dissipated by the unitary mass, p = \p\, L \s the integral scale, s is an integer exponent 
(typically, s = 2) and 

is the projector on the transverse degrees of freedom. 

Following the Martin-Siggia-Rose formalism [1] we introduce the Navier-Stokes density of Lagrangian 



C{v,w,P,Q,f)=w''{t,x) 



I - v"{t, X) + v^{t, x)-^v^t, X) 



'-Qit,x)^v'^{t,x) 



(8) 



where the field w" is the conjugated variable to the velocity field v" and the field Q is the Lagrangian multiplier 
related to constraint (4). The generating functional is given by the integral 

W (J, P)= j VvVwVQVf exp |i J dt d^x [C{v, w, P, Q, f) + J„t;"] 

-\j dtd^xd'yrF-^f^ (9) 

where J a are the components of the "external source" vector J. By integration over the statistical measure, 
and over the Lagrange multiplier Q, we obtain an expression which depends only on the transverse 
component Vt of the velocity field v. By decomposing the auxiliary field w in terms of its transverse {wt) and 
longitudinal (wl) components, w = wl + wt, the measure Vw factorizes into VwlVwt and the generating functional 
(9) reduces to 



Vwt'Dwl'Dvt exp <i dtd x 



d^xd^yw^Faffwi^ \ . 



(10) 
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Diagramatic strategies are usually applied at this level. We want to point out that one can go further by observing 
that also the transverse and longitudinal components of the auxiliary field w can be integrated out, yielding the 
equation 



W(J,P)= y VvTe-i'^''-^^+'^ ^^'^''^■^'''"^5 
where the action functional / has the form 



d 



1 dP 



(11) 



I{vt)= j dt(Px(fy 

Pa,'i\^-y\) 



d_ 

dt 



dt 



- vV'^j v^{t, x) + v^{t, x)dpv^{t, x] 
i^vA vl^it, y) + y)dxv^{t, y) 



(12) 



The computation of (11) would require to solve the constraint 



g d \ I dP ^ 



(13) 



In principle, this is a very difiicult task due to the nonlinear character ot the constraint. 

In the following section we show that we can identify a particular extremal solution, vt, of the functional (12). 
This solution is found to be independent of the stochastic source and, moreover, it satisfies constraint (13) for any 
constant value of the pressure. Accordingly, I{vt) can be interpreted as a large deviation functional (see eq.(l)) and 
the statistical nonequilibrium measure of the rsNSE can be effectively evaluated by integrating over the fluctuations 
around this extremal solution. It is worth observing that the entropy is related to the functional I{vt) by the relation 
[9] 



Sivr) = -^inf /(vt) , 

Z Vt 



(14) 



where the minimum is taken over all trajectories connecting vt to Vt ■ 

In what follows we are going to show that a suitable perturbative strategy can be applied for obtaining explicit analytic 
calculations of the statistical properties of the rsNSE. 



III. A QUASI-STEADY SOLUTION AND ITS STABILITY 

Any analytic approach aiming at the estimation of the generating functional (11) demands the identification of an 
explicit solution of the action functional (12). In practice, this amounts to solve the stationarity condition 
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Slivx) 



1 + ^V2 ) +d„v^{t,x) 



d_ 
di 



= . (15) 



+v^it,ff)dxvl^{t,y) 

We want to observe that for any arbitrary scalar field ^{t,x) a solution of the equation 

- i^V^) v^{t,x)+v^{t,x)dxv^{t,x) = d'^^t,x) , (16) 

is also a solution of (15). Since F^^Qx — y\) contains a projector on the transverse degrees of freedom we can fix, 
without prejudice of generality, the condition = 0. It is worth pointing out that, for what concerns eq.(16), this 
condition implies also that the longitudinal component of the nonlinear term vanishes, i.e. the solution Vj,{t,x) has 
to satisfy the additional condition 

df3(v^{t,x)dxv^{t,x)^ = . (17) 

Several different solutions can be found: among them, the only one unaffected by divergences in space and time is 
the following: 



v^{t, x) = — <l + e sin ( — ^ (b A a) • y ) > , with t>0 . (18) 



The U" are the components of the vector of velocity amplitude U {U = \U\) ,a=jj is the corresponding unit vector 
and b identifies a rotation axis. Both vectors U and b can be fixed arbitrarily. We assume also that the length-scale 
L is the same as the forcing integral scale defined in (7). This implies that solution (18) decays exponentially in time 
to the constant ^ with the rate td = 4L^/i/, which is the diffusion time scale. The dependence of solution (18) on 
the Reynolds number TZ can be made explicit by the relation 7?. = so that td = ^yV? /U'^. Notice that condition 
(17) is trivially satisfied by solution (18), because 

v^{t,x)dpv^{t,x) =Q . (19) 

Accordinlgy, Vj.{t,x) is also a solution of the diffusion equation (9( — i/V^) Vj.{t,x) = 0. There are two main con- 
sequences to be pointed out: i) as a solution of the linear diffusion equation v^{t,x) is unique, which is a crucial 
requirement for the large deviation approach; ii) the solution has to be defined only for positive times. 
Moreover, due to condition (19), the constraint (13) is trivially solved by P = constant. 
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For |;r| <^ L solution (18) approximates a linear shear flow: this is well known to produce instabilities for sufficiently 
large Reynolds numbers TZ. In this perspective, it is worth analyzing the dynamical stability of (18). For this aim we 
consider the perturbed velocity vector, whose components are: 

v" {t,x)= {t,x)+ Sv^ {t, x) (20) 

The perturbation vector 5v^ is assumed to be much smaller than with respect to any proper functional measure 
/i , i.e. \Sv^{t,x)\fi « \v°'{t,x)\fi yt and Vx . One can substitute (20) into (16) with d'^^ = 0, while assuming 
that it satisfies constraint (17). In the linear approximation one obtains an equation for Sv^{t,x), which can be 
solved explicitly by performing an expansion in the inverse Reynolds number TZ~^. As shown in Appendix A, the 
perturbation field vanishes and, accordingly, (18) is stable for sufficiently large times and Reynolds numbers and 
provided the following inequality holds: 

This inequality implies that for increasing values of TZ the band of unstable modes becomes thinner and thinner. 
As a consequence, solving the stability problem by expanding the solution of the linearized dynamics (Al) in powers 
of TZ~^ is consistent with this finding. Since condition (21) has been derived by assuming TZ large, it is not in 
contradiction with the Landau scenario for the origin of turbulence. 

In summary, v^{t,x) exhibits all the expected features of a physically relevant solution, which corresponds to 
stationarity conditions for the large-deviation functional. Accordingly, it can be effectively used for computing 
statistical non-equilibrium fluctuations of the rsNSE. In the next section we will exploit a saddle point strategy for 
performing explicit calculations from the generating functional. 

IV. PERTURBATIVE ANALYSIS OF THE GENERATING FUNCTIONAL 

All statistical properties concerning the rsNSE are contained in the structure functions which can be obtained 
by performing derivatives of the generating functional (11) with respect to the current J". An explicit calculation 
is unfeasible due to the nonlinear character of the action functional I{vt)- Since in the previous section we have 
identified the solution u^, we can tackle the problem by introducing the velocity field uf. = v^ — v^, which represents 
fluctuations with respect to v^, and by applying a saddle-point strategy. 
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Due to the translational invariance of the functional measure, the generating functional (11) can be rewritten as 

W{J) = j VuT e-5^("-r)+*/*'^'^''°"?. (22) 

A linearized expression for the action functional can be obtained by assuming that higher order terms in generated 
by the saddle-point expansion around the solution Vj. are negligible with respect to the functional measure Dut'- 

I{ut)= dt j SxSy [{dt - vVl)u^{x) + u^(x)apM?(i) + uf^{x)dpV^{x)\ F-^{\x-y\) x 

{dt - lyW^u^iy) + v^{y)dxu^{y) + u^{y)dx4iy)] ■ (23) 

We have also introduced the shorthand notation x = {t,x). 

Consistently with this perturbative approach, we can also assume that, at leading order, constraint (13) is still trivially 
solved by (18), i.e. the pressure P is a constant. 

In this way the action functional (23) has a bilinear form in the field u^, with coefficients depending on v^. In order 
to perform explicit Gaussian integration of the generating functional one has first to understand how the technical 
difficulties inherent such a dependence can be circumvented. The first problem that we have to face with is that, since 
(18) is defined only for t > 0, also (23) is defined for positive times. As we discuss in Appendix B, a standard procedure 
allows one to get rid of any singularity of the action integral that might emerge for t — > 0+. This is a consequence of 
the stucture of the linearized hydrodynamic operator appearing in (23). The second problem concerns the possibility 
of obtaining an analytic expression for the generating functional. To this aim one can exploit a perturbative expansion 
of (18) in powers of the inverse Reynolds number Actually, it is worth rewriting the solution (18) making explicit 

its dependence on the Reynolds number: 



v^(t,x) = — <l + e~^^'sm\ A T> \ } ,with t>0. (24) 

Using TZ~^ as a perturbative parameter, one can expand at all orders in TZ~^. When this expansion is substituted 
into (23) at leading order the action functional, in Fourier transformed variables, takes the form 

^^""^^ = / (04<(-^)^/(-P)^a/3'(P)<(^)4(P) + O . (25) 

We denote with u^(p) the Fourier transform of the field u^(x) with p= {po,p),Po and p being the Fourier-conjugated 
variables of t and x, respectively. We introduce the representation of the action functional in terms of the Fourier- 
transformed variables because this makes more transparent the diagonalization procedure required to arrive at the 
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final result. 

The hydrodynamic evolution term M^^{p)uj^{p) is given by the expression 

/' 1 \ r (bAuY 



r< (bAu) 
4 AvU 



> u'^ip), 



(26) 



where C = 



The next step in this calculation requires the diagonalization of the matrix Mp"(— "'^(p)M'^(p). Since by 
definition the factor [F'^l^{p)]~^ is proportional to the identity operator in the space of the transverse solutions^ we 
have just to diagonalize the matrix of the hydrodynamic operator M^^{p). 

The computation of the eigenvalues, Ai, A2 and A3 of M'^(p) deserves lengthy calculations sketched in Appendix C. 
Hereafter, we report the final form of the generating functional: 

y\){ri) = j J{H)'D4>T e~^/p<^T(-p)-F'"^(p)^pT(p)'^?'(p)+«/p'ira(-35)0T(p) ^ ^^2^) 

where we have used the shorthand notation Jp = J -(2t^ ^^'^ 

^ Xtip)Xrip) ^ 

Iplip) = A^(p)A2(p) , (28) 

^ xmMp)^ 

J{H) is the Jacobian of the basis transformation u — > J — > r/ engendered by the matrix H, which diagonalizes 

Mp°'{p). It is worth pointing out that the transformed vector (p^ip) still represents transverse components. Gaussian 
integration yields the following expression of the normalized functional in terms of the 77" source fields 



(29) 



In practice, the explicit computation of the structure functions can be accomplished by returning to the original 
representation, where the generating functional has the form 



>V(J) = e" 



-I If Ji{-p)F{p){HI-^H'^) mj^ip) 



(30) 



In the next section we are going to derive an explicit expression for the second-order structure function. 



iMore explicitly we have [F«''(p)]-i = P-'' {p)^^ {p) where F{p) = DqL^ {Lpf e'^'^"^ 
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V. SHORT-DISTANCE BEHAVIOR OF THE SECOND ORDER STRUCTURE FUNCTION 



The analj^ic expression obtained for the generating functional (30) allows one to obtain all the statistical information 
about the fluctuations around the basic solution Vj.. In this section we perform the explicit calculation of the second- 
order structure function of the velocity field u", defined as 

82= {\UT{t, f+ x) - UT{t, x)f) 

= {\{u'^{t,r + x)-u'i{t,x)){uTa{t,r + x)-UTa{t,x))\), (31) 

The brackets denote averages over the stochastic forcing. 

By assuming isotropy and homogeneity of the velocity field u", expression (31) is expected to assume the typical form 
of a scale invariant function 

52(r)=rf=F2(t,£) (32) 

Here r = |r| and L is the integral scale associated with the noise source. It is worth stressing that, at variance with 
fully developed turbulent regimes, here the assumption of isotropy and homogeneity have to be taken as a plausible 
hypothesis allowing for analytic computations. 



We want to point out that any exponent Cn must be independent of the basis chosen for representing the functional 
W. For the sake of simplicity, it is worth using (29) rather than (30) to obtain: 

6 5 \ f 5 5 



^^^^^ \i6r]^{t,x + r) iSr]^{t,x) ) \iSr]Ta{t,x + r) iSr]Ta{t,x'' ' ^^^^ 
As shown in Appendix D, it turns out that 52 (r) can be rewritten as follows: 

1 



7;=0 



where 



and 



S2{r) = --{hir)+h{r)). 



n=0 



r (2n + 4) 



Hr)= D,L^^-^ [f ^^Lpye-i-^r £^ ^x (e-™ - l) 



V 



E 

;=i,2 



C/2 

{1-x^y 



(l-a;2)3 



XS 



Ili^i (ELo i^ik - Sik) Fk {x, s, s)) 



y 



(33) 



(34) 



(35) 



(36) 
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The coefficients and the functions Fi, together with their arguments, are specified in Appendix D . 

The main contribution of the stochastic measure p^+*e~(^^'^rfp to the first integral in (36) comes from a narrow region 

of wavenumbers close to p, where the function p2+Sg-(Lp)^ ^las its maximum, i.e. 



Accordingly, the function ^^s^ p^ contributes to the integral by taking values close to ^^^^^ . 

Moreover, for p — p the sufficient condition (21) for the stability of small perturbations determines an upper bound 
for the Reynolds number: 

7^ < 4(s + 2) , (38) 

This implies that for sufficiently small TZ the wavenumber p is stable. Under this condition, the leading contribution 
in (36), consistently with the expansion in TZ^^, can be obtained by performing an expansion in powers of g^^^a ■ 
One finally obtains the complete expression of the structure function (see Appendix D for details) 

52(r)=-i(Ji(r)+/2(r)) 



(27r)2j/ 



^" .2vf, 1^„+ip/£+2n + 3\ 1 + S 2fs 2n + 4 /rx 
V 2 J r(2n + 4) Sf r(2n + 6) \lJ 



n=0 

for 1 < 7^ < 4(2 + s) , (39) 

At leading order in the distance r this expression is dominated by a dissipative contribution. 

We conjecture that this analysis can be extended to the parameter region defined by the condition TZ > 4(2 + s), 
where the statistically relevant wavenumbers can be unstable. As shown in Appendix D, in this case hir) has 
two contributions: one is again dissipative, while there is another one yielding the nontrivial scaling behavior r^/^. 
Specifically, the expression of 5*2 (r) for TZ > 4(2 + s) is found to be 

' Tru\ 4n ' r(2n + 4) \l) 

This expression is dominated by the term r^/^ for sufficiently small distances. Indeed, the crossover scale between the 
and the terms occurs at 

L^FTZ-i. (41) 

1j 
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In Appendix D we evaluate the constant F 0.6 and we report the expression of the numerical coefficient Co(S). The 
general expression of the coefficients C„(S) appearing in (40) has been omitted, because it has no practical interest 
for explicit calculations. 

It is a remarkable fact that S2 can exhibit the scaling behavior predicted by the K41 theory, which is assumed to 
hold (apart from intermittency corrections) when the velocity fluctuations are turbulent in the so-called inertial range 
of scales. This suggests that hydrodynamic fluctuations in a system at the very initial stage of instability development 
already contain some properties attributed to the developed turbulence regime. 

VI. CONCLUSIONS 

In this paper we have exploited the field-theoretic approach to reformulate the random forced Navier-Stokes problem 
in terms of the evaluation of a quadratic action. This has the formal structure of a large-deviation functional, 
describing thermal fluctuations of irreversible stationary processes. The crucial step for obtaining such a statistical 
representation is the integration over all longitudinal components of both velocity and the associated auxiliary fields. 
With respect to the standard formulation which yields usual diagramatic strategies, we perform one more field 
integration. The positive definite kernel, which connects the hydrodynamic evolution operator in the action functional, 
is the inverse of the forcing correlation function. 

In terms of the action functional, the knowledge of the whole velocity statistics reduces to the computation of functional 
integrals. However, due to the intrinsic nonlinear character of the hydrodynamic operator several technical difficulties 
have been solved for performing analytic calculations. In particular, one has to to introduce suitable approximations. 
In order to obtain an analytic expression of the generating functional we have identified a solution around which 
we have linearized the hydrodynamic evolution operator. We have also introduced a velocity field which represents 
fluctuations with respect to this solution. A perturbative expansion in the inverse Reynolds number finally yields the 
wanted result. 

In principle, from this analytic treatment one can obtain all relevant statistical information about the rsNSE by 
computing any velocity multipoint structure function. In this paper we report only the explicit calculation of the 
two point second order moment of the velocity field. As shown in the Appendices, the algebraic manipulations needed 
for obtaining the final result are far from trivial also in this simple case. 

In fact, in this paper we aim at understanding whether fiuctuations at the early stage of their development (accord- 
ingly, we dub them as pre-turbulent fiuctuations) already contain some important features of developed turbulence. 
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We are interested, in particular, to characterize the scale invariant properties of such fluctuations. In this respect, we 
find that they are organized at different scales in a self-similar way. Remarkably, the scaling exponent coincides with 
the dimensional prediction of the Kolmogorov 1941 theory [11], valid for developed turbulence regimes. Whether or 
not such exponent is a genuine reminiscence of the developed turbulence phenomenology needs further investigations. 
Unfortunately, the complexity of the derivation leading to the K41 scaling law does not allow us to identify precisely 
the very origin of such a dimensional prediction. We can however argue a relationship between the observed dimen- 
sional scaling and the conservation laws (for momentum and energy) associated with the two eigenvalues of the matrix 
appearing in the action functional (25). 

Finally, it is worth observing that the dimensional scaling law emerges for a particular choice we did for the pressure 
field: fluctautions have been restricted around a solution for which the pressure is constant. Unfortunately, owing to 
the fact that the analytical treatment is not duable in the general case, we cannot substantiate the fact on whether 
the dimensional prediction we found is not a consequance of our particular choice for the pressure fields. 
At least three scenaries might be possible. Firstly, pressure field does not affect neither the leading (dimensional) 
scaling law nor its prefactor. It only affects the subleading scaling contributions. In this case our simplification would 
capture the relevant physics of the problem. The second possibility is that the leading scaling law does not change 
but this is not for the prefactor. The last possibility is that pressure changes the (domensional) scaling law giving rise 
to intermittency corrections. Unfortunately, at the present stage of our knowledge, we are not in the position to select 
one scenary among the three we have pointed out. Further investigations are needed for this aim, which probably call 
to deep numerical investigations of the system under consideration. 

We want to conclude by outlining some open problems and perspectives. A first question concerns the physical 
relevance of the solution (18) around which we linerize the evolution operator. It represents a shear-like solution, 
which is a well-known generator of instability. Moreover, its unicity and stability properties seem to indicate that this 
solution can play a major role in the determination of stationary nonequilibrium fluctuation statistics to be attributed 
to the rsNSE. As a mathematical object, it exhibits all the wanted features that one would like to attribute to such a 
solution. On the other hand, the authors have not yet a physical intuition for its relevance and aim at making some 
future progress in this direction. 

Another interesting point to be tackled concerns the computation of the third-order momentum of the velocity cor- 
relators. In this case the predictions of our approach could be compared with the 4/5-law, which is one among the 
very few exact results of turbulence theories. 
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Finally, the extension of our results to other classes of transport problems, including passive scalar advection, could 
provide a better understanding of the basic mechanism at the origin of the observed scaling behaviors. 
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APPENDIX A 



In this Appendix we perform the stability analysis of the solution by the linearized equation 



d_ 

dt 



z^V^^ 5v^{t, x) + v^it, x)^-p6v'].{t, x) + 6v!^{t, x)-^-^v^{t, x) = 



(Al) 



with the constraint 



d 
dx-y 



(^v!^{t,x)-^dv'^{t,x) + 6va^{t,x)-^v'^{t,x)^ = 



In Section III we have already observed that is a quasi-steady solution for a time t <^ td = ^f^- The Fourier 
transform of eq.(Al) with respect to the space vector x yields: 



dt 



Sv^ (t, kj - lyk^Sv^ (t, k) + ^k- U6v^ (t, + \e~^ jc/'^fc, 



Sv"^ i, fc - C- 



-Sv^ {t,k + C 



b/\U 



+ U°'C- 



s4 t,k-C^-^]+s4 t,k + C^-^ 



= 0, 



(A2) 



where C = ^(~g)2 " performing a perturbative expansion up to second order in the parameter H. ^ , one obtains 
the system of equations 

d 



k ■ USvl 



'T(O) 



(t,k) 



§-^Sv^ii) (t, k) + vkHv^^^^ (t, k) + '-k- U6v^^,) (t, , 

1 (b/\u) 
(t,k)--U"C^ 



1. (^^^) 

-k-UC 



4^dk~0^^°^ 



^5vt.. {t,k) 



(A3) 



(A4) 



T(2) (^1 ^) 



'T(2) ^) ~'~ 2^ ' ^ ^^T{2) (^5 ^) 
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^ (bAu) ^ (bAu) 

tk-uc^—f^^Sv^,^. {t,k) - ic/"C^-— (t,k] 
2 4:vH dkf) \ J 2 4:vn ^(1) V J 



This system of equations yields the perturbative solution 



6v^ 



(t, k) = e-i^'^'+i^-^^y I F("o) (k) + (k) 



k-U 



(bAU)- Vfei^("o) ^(0) (k) t 



SuTZ 

(bAliy kF^o^ (kj vi 



where the functions F's are determined by the initial conditions: they are found to be of 0(1) for any k. 



(A5) 



(A6) 



The exponential term in front of (A6) makes the perturbative solution vanish in the limit of large time provided 
the perturbative series contained in the curly brackets does not diverge faster in such a limit. This requirement can 
be translated into the following spectral condition 



C/2 



-fc^ > 1 . 



(A7) 



This inequality indicates that the instability of solution (18) may originate only from sufficiently small values of the 
wave-number k. 



APPENDIX B 

As shown in Section III the solution of the hydrodynamic operator in the action functional (12) is defined for 
t> Q. Accordingly, it breaks Galilean invariance, thus giving rise to the well-known Doppler effect, i.e. k^ — > ko+^k-U. 



Moreover, since in Section IV we evaluate the action functional by applying a saddle- point expansion around w^, 
the approximated expression (23) contains a time integral that has to be restricted to t > only. This amounts to 
assume that the action should be identically zero for f < 0. Accordingly, one cannot exclude the possibility that a 
singularity in the time integral may originate at f = 0. 

In this appendix we want to show that one can easily exclude the presence of any singularity by passing to a 
Fourier-transformed representation of the action functional (23): according to a standard field-theoretic technique 
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the addition of a small immaginary part to the frequency appearing in the Fourier-transformed integral allows one to 
control its regular behavior for t — > 0+. 

For the sake of clarity, we present this procedure only for two of the terms appearing in (23). Actually, one can 
easily realize that the procedure can be extended to all the terms: we just report the final result, thus avoiding the 
writing of lengthy formulae. 

Let us consider the term 



h = 



£ dt J d?x J d'^y^uUt,x)F-'"0{\x-y\)^^u^{t,y) (Bl) 



In principle, the integral in the time domain is ill-defined. We can pass to Fourier-transformed variables and rewrite 
it as follows: 

The time integral can be regularized by adding a small immaginary part ie to the frequency component and the 
integral Ii is transformed into 

- / (0 /_r ^ /r s r ^^^^^^^^^^^^^^^^^ -^^) 

d^k f+^dko f+^dq„ k,q, . ^ ^ 

Urp[ko, k)Uj.[qo, -k) 



(2^)3 7_ 2nJ_^ 2.z(fco + 9o + ie)™'^^"^ 

d^k f+°^dk„ /•+°°dgofco(<?o-fco)- 



hpr^ ^/ ^5 M^(fco,fc <(go -fco,-fc) . B3 

J {2TTy 27r 2tt i{qo+ie) 



By performing the limit e — > 0+ one obtains 



-^1= ~ / 7^~T3* / ^ Pf ^—ko{qo-ko)u^{ko,k)u^{qo-ko,-k) 

-^S{qo)ko{qo - /,;o)?iT(/eo- ^)mt(9o - fco, -k) 
-oo 27r 

= T^T^ / ^ / u^{ko,k)ut^{qo- ko,-k) 

J y^'^r J^oo 27r [27r qo 



+ ^^o^t(^o, k)u^{-ko, -k) 



(B4) 



In this equation P denotes the principal value. The nontrivial part to be computed is contained in the square brackets. 
One has to consider that the fluctuations Uj.{t,x) become negligible for scales smaller than the Kolmogorov scale. 
Since they are defined for t > and the time integral is singular in t = 0, we have that its Fourier-transformed 
representation should exhibit a unique singularity at infinity, where it vanishes for Im qo < 0. One can write: 
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27r 



-P 



ctgo ^ u^{ko,k)utp{qo - ko,-k) 



90 



k^u^{ko,k) 

2lT 



+ 00 -/S 



dqo 



u'j,{qo - ko, -k) 



-kQu'^{ko,k)uij,{-ko,-k) . 



Making use of this result, one can easily conclude that (B2) can be written as follows: 

dkod'^k 



(B5) 



k'^u^{ko,k)F-'^"^{k)u^{-ko,-k) • 



(B6) 



Now, let us consider one of the terms of (23) which exhibits the Doppler effect in its Fourier-transformed representa- 
tion: 

/f Q 
d^x J d^y-u^{t,x)F-''"^{\x - y\)v^{t,y)dx4{t,y) 



c+00 

+ I dte' 



■I 

Jo 



kx + C 



bAU 



+ I - ^ J I 50, -k + G- 



-2j 2^y_< 



AvTZ 

i 



d^k f+°° dko f^^dqol 

2tt 2n \ko + qo+ 



F-^l^{k) 



1 



(fco + 90 +«4^) 

We expand the solution up to first order in powers of TZ~^ and wc obtain the final expression: 



(B7) 



1. = -, 



1 fdkofk_^ r,M^-u) [~/3 



(27r)4""^(^°'^)>"/3(A:) 



Urp{—ko, —k) + iC 



.r<{bAU)x d 



M^(-fco,-fc) + 0( ^ 



As in the previous case, one can regularize the integral in t = by performing the limit e — > 0+. By applying this 
procedure to all of the remaining terms in (23) one arrives at the final expression (25). 



APPENDIX C 

In this Appendix we sketch the calculation of the eigenvalues of the matrix M^{p) defined in (26). In fact, the 
perturbative expansion of the solution (18) in powers of^ induces an analogous expansion for this matrix. Formally, 
one can write 

M = M(o)+M(i) + ... (CI) 
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where 



i(po + ^P- U] +i^p^ 



bAU 



Mf,) -6^,jp. U^^d,^ - jU-^^ . (C2) 

The matrix Af^f (p) acts on the two-dimensional space of the transverse functions and on the one-dimensional space 
of the longitudinal functions. Only the transverse degrees of freedom are physically relevant. 



A complete orthonormal basis in is given by the vectors 

(bApY 



gip) {bApy - f{p) (uAp 
^f{pWf{p)Kp)-9^{p) ' 



n3«= ^ , (C3) 

where 

f{p) = hV-{b-p)\ g{p) = {b-U}p' -{b-p){U-p), h{p) = uV-{U-P)^ . (C4) 

n" and span the transverse subspace, while Ilf spans the longitudinal one. In analogy with (CI), also the 
eigenvalues of M^{p) can be represented by a perturbative expansion in powers of^, namely as 

A° = A^o) + ^(1) + - where a = 1, 2, 3 . (C5) 

The zero-order eigenvalues A°q^ are degenerate and have the form 

A^o) = (^(po + ^p-C/)+V) . (C6) 

The evaluation of the first order corrections A^j^^ requires the diagonalization of the matrix with elements = 
(nj,M(i)nj), = 1,2,3). After some simple but lengthy calculations one finds 

A[i)= ^ (^M(i)n + M(i)22 - \/ (M(i)ii + M(i)22)' + 4.M^i)2iM^i)i2^ , 
>\i)= \ + M(i)22 + \] (M(i)n + M(i)22)^ + 4M(i)2iM(i)i2^ , 

A^i)= M(i)33 , (C7) 

with 
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(O (bAU)-p 



C ((bAU)-p)(b-p)g{p) 



M, 

C {{bAu)-p)(h-p 



M(i)i2 



16z/7e/(p)//(p)/i(p)-g2(p) 



6- 



C/)ff(p) + 2(p-C/)ti;(p)-C/V(p) 



fi ((bAU)-p) ^\ r/ \ / \ 



Mn 191= , 

M(,,33=-^fc^i(pA[7).,^) , (C8) 

where we have introduced the further definition: 

w{p) = b'^{p-U)-{b-p){b-U) . (C9) 

Without prejudice of generaUty, we can specify the geometrical structure of the flow. For the sake of simphcity, we 
assume that the vector r (i.e. the Fourier-conjugated variable of p) corresponds to the polar axis and that the vector 
b is orthogonal to both r and U. With this assumption the two physically relevant first-order corrections to the 
eigenvalues are 



f/2 



A(i)= iQi^-ji {^^^^u cosOu [cos^ (j)u + cos {2{(j)u - (j)))] sin^ 6 

+ cos^ 9u sm29 cos{(j)u -<l))} ■ (CIO) 

Since A^^^ is associated to the longitudinal part, it does not play any role in our calculations. 

APPENDIX D 

In this appendix we aim at reporting the main calculations needed for obtaining an explicit expression for (33). 
According to the perturbative approach discussed in detail in Appendix C, S2{r) can be written as follows: 

W ~ -2 / - 1) ± — , . (Dl) 

The eigenvalues A"^^ which appear in this equation have been computed up to first order of the perturbative expansion 
in TZ~^. Notice that the sum is restricted to the first two eigenvalues {a = 1,2), which correspond to the transverse 
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components of the velocity field. Actually, the third eigenvalue, corresponding to the longitudinal components of the 
velocity field, is ineffective for our calculations. 
Explicit integration over po yields 

With the particular choice performed in Appendix C for the geometrical structure of the flow, S2{r) can be expressed 
as the sum of two terms: the first one is associated with the null eigenvalue Aj^^, while the second one depends on 
the nonzero eigenvalue Namely, 

S^{r)=-l{h{r)+h{r)) (D3) 

By considering the explicit expressions of the statistical function F{p) and of the eigenvalues A"^^ (see eq.(ClO) ), one 

has 

I (r)= DL^ I ^ {LpYe-i^v)' 

'^^ ' ° J i2nfp2 + j^^[2smeuCOseusm^9cos^(j) + cos^eusm2ecos(j)] 

(D5) 

In the r.h.s. of this equation we have also exploited translational invariance for applying the transformation {(j)ij — <p) — > 
—(j). The analytic calculation of (D4) is obtained by a standard procedure: 



27r2 ^ (2n)!(2n+ 1) 

= f (_i)n+i£(flXl!l (D6) 

(27r)^ r(2n + 4) \l) ^""'^ 

For what concerns we first perform the integration over the variable 0, namely: 

h{T) = DoL^ ^(Lp)^e-(^f)' d{cose) {e'^^-°^' - l) /q (D7) 



where 



Jo 



1t! 



P^ + 1^3^ ""^^ ^ '^^^ ^ '^'-'^^ 2^ '^'^^ ^] 

32i/27e r zdz 

az^ + bz^ + cz'^ + bz + a ' ^ ' 
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with z = e^'^ and the integration is on the unit circle 7. The coefficients a,b,c are given by 

a = sin 9u cos 9u sin^ 6 , 6 = 2 cos^ 9u sin 9 cos ^ 
32^/2-^ 



c = 



(72 



-p + 2 sin 9u cos sin ^ . 



(D9) 



The evaluation of the integral (D8) requires the knowledge of the root of a fourth-order algebrical equation. By 
exploiting the Euler method [14] we end up with the expression 



Zi — Zi I X, 



X3 



(1-X2)S 



^ / 8l/27^ 2 ^ A iv. 



.1=0 



X3 



(l-a;2)3 



i = 1,2,3,4 . 



The following definition has been adopted: 



with 







m 


Ef 


"81 4 


12 





(1 - ,X2) 

64 1 - a;2 ^v^U 



81 . .x2 
2 ^ 2 (l-x2) 



90 



E2 



C/2 



-r 189— 



382 



- 120 



8:^27^ \ ( 504 
C/2 ^ j 

X I e' [1 + (1 - 4 X 27 /i) 

8l/27^ n 



S (l-x2)2 5(1 -.r2, 

47 E2 \ /8i^27^ 



S2 (1^x2)2 S2a;2(l-a;2) 



^ 2:2 

3 



C/2 

1 - (1-4 X 27 



32 E2 



S3 a;2 (1 - a;2)2 

4 4 
1 S3X3 



(1-^2) 



2^ 6 



3 E3a;3 



f/2 



-P 



3E3Sx3(l-a;2)3 



(DIO) 



a; = cos 6 , S = sin cos , 5] = cot ^[/ 
/ . . , \ 



1 


1 


1 


1 


-1 


-1 


-1 


1 


-1 


-1 


-1 


-1 



Here e is the cubic root of unit: e 



16 + 30 S2 



111, 



. The explicit expression of the function h follows: 



l-a;2 



(l-a;2)2 
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8t/27^ p2 i7E2a;2 + 3(1 - a;2 
i72 5 (l-a;2)2 



(Dll) 



(D12) 
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+48 



U^E J (l-x2 



H 3 



-128 + 81 S' 



. x' . 81^6 

(1-X2)2 ^ 2 (1 - ,t2)3 



-180S 



2 , 8^.27^p2 / _ s4,t4 



+ 



l-a;2 
8^/27^p' 



+ 



2^2 



378- 



[72 S V"'^(l-x2)-' 

£2^2 „ . 1 



764 



240- 



1 



1008- 



94- 



(l-x2)2 l-a;2 

8^/27^p2^3 ^ 



64 



[/2 s; V (l-x2)3 (l-a;2)2; \ E J {1 - x^)^ 

Only the roots zi and ^2 are included into the unit circle, therefore (D7) becomes 

^ , , ^ 3 32l.27^ p2^p 
J2(r)= DoL" 



E 



(D13) 



f/2 7o ^,,),iFpy^-''^' Jj^i^---^) 



l-x 



2\3 



ELo {x, ^p'; S, h) + iE-^ 



(l-x2)3 



i=l,2 



a;3 



(D14) 



U^^l (ELo (^^fc - {X, ^p2. S, S)) 

As we have already observed in Section V, only the values of the variable p around p = give a significant 

contribution to the integral in (D14). We observe that ^jj^p^ — > ^^^^^^ and the stability condition (21) imposes: 



l<7e<4(s + 2) . 

The evaluation of the leading terms is then possible by performing an expansion in the parameter -g^r^p '^ 
that, by virtue of (D15), is smaller than unit if C < '\J ^2^- 

For C > V the contribution to the integral rapidly vanishes. For 1 < 7?. 4; 4(s + 2) we obtain 



(D15) 

8 



S,{r)=--ih{r)+h{r)) 



2 

r 



n=0 



s + 2n + 3 



1 + S 



2-T5 2n + 4 



r(2n + 4) eI r(2n + 6) 



(D16) 



By extending the validity of our calculations to 7^ > 4(s + 2), wc have ^jj^p'^ < 1 for ^ < y As in the 

previous case, we expand (D14) in power of the parameter ^("^ < 1 and we obtain: 



/2(r)- EDoL' 



+^E 



(27r)2 
(l-x2)^ 



1 r 1 + 8.^2 , 



Em=o s/m-Fm (a;, 0; S, S) 



" i^[^2 Vn^^/ (ELo - ^ik) Fk {x, 0; S, S)) 



H — C 



d 


E1=o simFm {x,y;T,,^) 




Ui^i (ELo i^ik - s^k) Fk (x, y; S, 



y=o 



(D17) 



Two different terms, /2^(r) + /^(r) = hir), can be identified in (D17). The evaluation of the first term is straight- 
forward: 
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n+1 



ED, - (-1) 
2(27r)2 A.r(2n + 4) 



s + 2 + 2n 



s + 5 + 2n 



(i) 



2n 



7e 



(D18) 



The evaluation of the second term is more cumbersome. The leading term can be recasted in the form: 



(D19) 



The coefficients Ai are S-dependent numerical constants. The first two of them are given by the expressions 

Ao(E)= ^ 



16\/3(l-sinf)cos(itan-i\/26) ' 

sm ( 1 

Ai(S)=-- 



65sin(|tan-i V26) 



(D20) 



512^26 cos2 (i tan-i ^26) ' 

The exact form of these coefficients is however irrelevant for our analysis. Some tedious standard calculations yield: 



\o/ n=0 ^ 

where the coefficients C„(E) depend on the constants Ai. For n = one has 



(D21) 



(D22) 



The comparison between /;f (r) and /^(r) indicate that a crossover between the corresponding scaling behaviors occurs 
at 

-0.0336- 0.1127 cot^ Ou 



f 



2 X 8.328v^- 



(D23) 



2 + sin 9u cos 6u 

For the perturbative expansion in ^ to be meaningful, the parameter 9u must have a value close to |-. This implies: 



r ~ FTT^L, with F ~ 0.6 . 
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